Adsorption of charged particles on an oppositely charged surface: Oscillating 

inversion of charge 
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Adsorption of multivalent counterions on the charged surface of a macroion is known to lead to 
inversion of the macroion charge due to the strong lateral correlations of counterions. We consider 
a nontrivial role of the excluded volume of counterions on this effect. It is shown analytically 
that when the bare charge of macroion increases, its net charge including the adsorbed counterions 
oscillates with the number of their layers. Charge inversion vanishes every time the top layer of 
counterions is completely full and becomes incompressible. These oscillations of charge inversion 
are confirmed by Monte-Carlo simulations. Another version of this phenomenon is studied for a 
metallic electrode screened by multivalent counterions when potential of the electrode is controlled 
instead of its charge. In this case, oscillations of the compressibility and charge inversion lead 
to oscillations of capacitance of this electrode with the number of adsorbed layers of multivalent 
counterions. 
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PACS numbers: 82.70.Dd, 87.16. Dg,87.14.Gg 



I. INTRODUCTION 

Adsorption of charged particles on the surface of an op- 
positely charged macroion is an important problem with 
broad applications in many areas of science. In a wa- 
ter solution, double helix DNA, actin, charged colloids, 
charged membranes or any charged interfaces can play 
the role of the macroion. Charged particles, or coun- 
terions, in solution can be ions, small colloids, charged 
micelles, short or long polyelectrolytes. Mean field the- 
ories based on the Poisson-Boltzmann equation and its 
linearized version, the Debye-Hiickel equation, have been 
used to study such screening problem. There are, how- 
ever, several new phenomena in solutions containing mul- 
tivalent counterions which cannot be explained using 
standard mean- field theories. The most notorious of all 
is probably the charge inversion, a counter intuitive phe- 
nomenon in which a macroion strongly binds so many 
counterions that its net charge changes sign. This can 
be thought of, theoretically, as overscreening. It cannot 
be explained using Poisson-Boltzmann theory because, 
in a mean field theory, screening compensates at any fi- 
nite distance only a part of the macroion charge. Charge 
inversion recentlirlras aiJpacted significant attention.ETO 

It was shownme-BBthat charge inversion is driven 
by the counterion correlations which are ignored in mean 
field theoriesE2l. When multivalent counterions condense 
on the surface of a macroion to screen its charge, due to 
their strong lateral repulsion they form a two-dimensional 
strongly correlated liquid. To see when correlation is im- 
portant in this liquid, one first defines a dimensionless 
parameter F which measures the strength of the interac- 
tion energy between counterions in units of the thermal 
energy: T = {Zefi-nnf^/DkBT, where D ~ 80 is the 
dielectric constant of water, Ze is the charge of a coun- 
terion (for convenience, we call it a Z-ion) and n is the 
two-dimensional concentration of Z-ions condensed on 
the macroion surface. Without loss of generality, we as- 



sume through out this paper that Ze is positive and the 
macroion is negatively charged with the surface charge 
density —a. For a strongly charged macroion with sur- 
face charge density —a ^ —le/nm?, F ~ 1.2, 3.5, 6.4 
and 9.9 at Z = 1, 2, 3 and 4. Thus, for Z > 3, F is a 
large parameter of the theory and the two-dimensional 
system of Z-ions is a strongly correlated liquid (SCL). It 
has short range order very similar to a two-dimensional 
Wigner crystal (See Fig. (|l|)). At the same time, large 
parameter F also guarantees the layer of Z-ions at the 
macroion surface is effectively two-dimensional. Indeed, 
each Z-ion, in the uniform field 2'Ka of the macroion sur- 
face, moves within a distance DUbT /2-KaZe ~ 0.52yl/F 
from the macroion surface. At large F, this distance is 
much smaller than average distance 2^1 between Z-ions 
in the direction parallel to the macroion surface. This 
makes the Z-ion layer effectively two-dimensional. 
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FIC. 1. The two-dimensional strongly correlated liquid 
(SCL) of Z-ions on the macroion surface with the short range 
order of a Wigner crystal. Black disks are Z-ions with radius 
a. The lattice constant of the crystal is 2A = y2/\/3n. 



The -c.orj'p.lation physics explains charge inversion as 
foUowsB^EJO. When a new Z-ion comes to the macroion 
which is already neutralized by the Z-ion layer, it pushes 
other ions aside and creates a negative background charge 
for itself. In other words, it creates an oppositely charged 
image (or a correlation hole) in the Z-ion layer (Fig. 
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ga). Due to the attraction from this image, the Z-ion 
sticks to the surface and overcharges it. The overcharg- 
ing saturates when the positive inverted net charge of 
the macroion is large enough to counterbalance the above 
mentioned attractive force. In a more quantitative lan- 
guage, correlation induced attraction to the Z-ion layer 
can be explained as a result of the negative chemical po- 
tential ^JLscL of SCL, whose magnitude is much larger 
than keT (|mscl| ~ TkeT > keT). 

In this paper, we would like to investigate the problem 
of charge inversion further by considering the effect of 
the excluded volumes of Z-ionSf— This effect was studied 
in Poisson-Boltzmann approachEll. Here we show that in 
the correlation theory, this effect can lead to qualitatively 
new consequences. When the charge of the macroion in- 
creases, A decreases and F increases, the correlation be- 
tween Z-ions becomes stronger. In other words, l/iscil 
increases and charge inversion is stronger. However, 
the excluded volume has a dramatic consequence when 
A ^ a: the pressure of the SCL liquid diverges and it 
becomes effectively incompressible. 



a) b) c) 

FIG. 2. The origin of oscillations of charge inversion, a) 
The adsorbed Z-ion layer is only partially filled. A new Z-ion 
pushes other ions aside to create a correlation hole for itself. 
The Z-ion is closer to the negative surface charge of this cor- 
relation hole than to the positive neighbouring Z-ions and 
therefore is attracted to the surface, b) The layer is filled and 
correlation hole disappears, c) More than one layer is filled. 
The correlation hole of the approaching Z-ion is created in 
the top layer only. 

In this case, the correlation image and charge inversion 
disappears (Fig. ^b). If the macroion charge increases 
further, however, a second layer of Z-ions is created and 
charge inversion recovers (Fig. |^c). Here everything is 
repeated again but with the role of a played by the net 
charge of the macroion and the first full layer. Conse- 
quently, if one plots the net charge of the macroion ver- 
sus its bare charge, one sees an oscillation of charge in- 
version: every time a new layer of Z-ions forms on the 
macroion surface, charge inversion appears and steadily 
increases with increasing macroion charge then drops 
quite abruptly when the layer starts to become full (see 
Fig. H). As a result of these oscillations, the maximum 
value of the inverted net charge is limited, no matter 



how large the charge of the macroion is. To our knowl- 
edge, this is the first time this periodical behaviour is 
suggested. 
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FIG. 3. Schematic plot of the macroion net surface charge 
density a* as a function of the absolute value of its bare charge 
density a. Every time a layer of Z-ions is almost full, charge 
inversion disappears and then recovers when the next Z-ion 
layer builds up. The dashed line which is described by Eq. 
(H) is the maximum charge inversion one obtains if the effect 
of excluded volume is neglected [a = Q). 

This paper is organized as follows. In Sec. II, we de- 
rive an analytical formula for the chemical potential of 
the Z-ions condensed on the macroion surface taking into 
account the excluded volume effect. This chemical poten- 
tial is then used to calculate the inverted charge of the 
macroion. In Sec. Ill, a Monte Carlo simulation of a 
small system is used to verify the analytical predictions 
of Sec. II. In Sec. IV, we consider the role of charge 
inversion in the case of screening of a metallic electrode, 
whose potential is tuned instead of surface charge density 
a. We show that, in this case, correlations lead to the 
oscillations of the macroion capacitance. 



II. CHEMICAL POTENTIAL OF Z-IONS WITH 
HARD-CORE REPULSION. 

To understand the electrostatics of the system, let us 
start with a two-dimensional system of point charges 
(a = 0) with concentration n on a neutralizing back- 
ground charge. Conventionally, this system is referred to 
as a one component plasma. Excluded volume effects will 
be considered later. At large parameter F, the plasma is 
a strongly correlated liquid (SCL) whose short range or- 
der is very similar to a Wigner crystal (WC). Specifically, 
the energy per Z-ion escLin) of a SCL can be very well 
approximated by that of a WC, which in turn can be 
approximated by the energy of a Wigner-Seitz cell be- 
cause quadrupole-quadrupole interactions between neu- 
tral Wigner-Seitz cells are very small. This gives 
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A more accurate calculationE3 gives ewcin) a slightly 
higher value: 



ewc{n) ~ -imZ^e^/AD 



(2) 



Knowing £wc{n), the chemical potential of the WC can 
be easily calculated: 
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Thus, nwcin) is negative and much larger than fcsT. It 
acts as an additional attraction to the macroion surface 
which in turn leads to the charge inversion effect. The 
temperature correction to nwc{n), or in other words, the 
difference between u^cLin) and ^iwcin), has been cal- 
culated numericall}^ for a wide range of the parameter 
r and indeed is very small and can be neglected. For 
example, it is equal to 11% at F = 5 and 5% at F = 15. 

Now let us return to our problem of finite size Z-ions. 
As the detail derivation in the appendix shows, the chem- 
ical potential of a Z-ion with hard-core repulsion can be 
written as the sum 
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where /xq is the chemical potential of a system of neutral 
hard discs of radius a at the same concentration n. In 
writing down Eq. (^, the electrostatic contribution to 
fi{n) has been approximated by fJ-wcii^), in exactly the 
same way fiscLin) was approximated by fiwc'in) above. 
The temperature correction to fiwcin), in this case, is 
even smaller than that for the SCL of point charges be- 
cause due to the hard-core repulsion, two Z-ions are never 
at a distance less than 2a from each other. Further- 
more, this correction decreases when A decreases making 
the approximation by fJ-wcii^) increasingly better and 
asymptotically exact when A — > a. 

To proceed further, one needs to calculate the chem- 
ical potential Ho{n). The hard disc and hard sphere 
systems have been extensively studied both analytically 
and numerically providing significant insights into solid- 
liquid phase transition and the melting process (See Ref. 

and references therein). In this paper, to emphasize 
the effects of excluded volume, we would like to repeat 
here simple calculations using the free volume approxi- 
mation which is asymptotically exact near the close pack- 
ing limit. 

To calculate /j,o, one notices that because of the hard- 
core repulsion, the area s available per each disc is smaller 
than the area 2y/3A'^ of a Wi gner-Seitz cell (See Fig. H). 




FIG. 4. The center of the central Z-ion is confined in the 
gray area by the hard core repulsion from its six nearest neigh- 
bours on the hexagonal lattice (the black dots). Near close 
packing, this area is significantly smaller than the area of the 
hexagon itself. 

A simple calculation shows that near the close packing 
limit, s = a{A-a)'^ < 2^/3^^, where a ~ 4.62. To calcu- 
late the free energy of the quasi two-dimensional layer of 
Z-ions, one needs to add a third out-of-surface dimension. 
The typical distance x of the out-of-surface motion of a 
Z-ion is the distance at which the energy cost for moving 
the Z-ion away from the plane is of the order of ksT. 
Solving 2'KaZex/D ~ ksT gives x = DkBT/2TTaZe = A 
which coincides with the Gouy-Chapman distance (note 
that this estimation of x is consistent because, for F ^ 1, 
K A, all other Z-ions are far apart and we can use 
the planar electric field 21X0 jD for the calculation of x.) 
Now, one can estimate the free energy of the hard disc 
system as the free energy of an ideal gas with concentra- 
tion 1/sA: 



Fo = -NksT ( In — + 1 



(5) 



where N is the number of adsorbed Z-ions and vq is the 
normalizing volume which can be thought of as the vol- 
ume of one water molecule. The exact value of wq is not 
important because it drops out of all final results. 

The chemical potential of adsorbed Z-ions is, corre- 
spondingly: 
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From Eq. (||), one can see that when A — a <C a, this 
chemical potential is positive and very large what negates 
any gain from the electrostatic part \iwc- Charge inver- 
sion therefore disappears in this limit and will not re- 
cover until the second layer builds up when the macroion 
charge grows. 

Using the correlation hole concept mentioned in the 
introduction, one can easily understand why charge in- 
version disappears. The second term in Eq. (|^), which 
is the dominant term when A approaches a, is nothing 
but po/"^, where po = —dF^/dS = nksTA/ {A — a) is the 
two-dimensional pressure of the hard disc system near the 
close packing limit, S is the total area of the macroion 
surface. The divergence of the pressure leads to the zero 
compressibility of the Z-ion layer when A ^ a. When 
a new Z-ion approaches this layer, it cannot put other 
Z-ions aside to create its correlation hole to attract to it 
(Fig. ||b). Therefore, charge inversion disappears. 

It should be noted that, even though Eq. (^ is ex- 
pected to be asymptotically correct only near the close 
packing limit {A — > a), where the hard disc system is 
known to be in the solid phase, it is actually quite good 
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as an approximation in the whole range of A. Indeed, 
Monte-Carlo simulations and molecular dynamics calcu- 
lations (See Ref. ^5|) show that the pressure obtained 
using the free volume approximation gives an excellent 
result for the solid phase of the hard disc system and a 
good result for the liquid phase. The biggest discrep- 
ancy occurs around the liquid-solid phase transition is 
only about 20%. Combining with the fact that, for large 
r, the electrostatic contribution fxwc is the dominant 
term of fi everywhere except near close packing, we can 
reliably use Eq. for the Z-ion chemical potential. 

Let us now calculate the net charge of the macroion 
which is the total charge of the macroion and the Z-ion 
layer on its surface. To do so, one has to balance the 
electro-chemical potential of the Z-ions at the macroion 
surface with the electro-chemical potential /it, of Z-ions in 
the bulk solution. If we assume the concentration of the 
Z-ions in the bulk solution is c, this balance in electro- 
chemical potential reads: 

li + Ze^{0) ^ fib — kBTln{cvo), (7) 

where ip{0) — 2'Ka*rs is the averaged macroscopic elec- 
trostatic potential at the macroion surface due to the net 
charge density cr* — —a + Ze/2^/3A'^ of the macroion 
surface with Z-ions condensed on it, Ts is the screening 
length of the solution. On the right side of Eq. (^, the 
bulk solution (where the electrostatic potential vanishes) 
is assumed to be ideal so that /if, is approximated by the 
ideal gas chemical potential kBTlTi{cva). 

Substituting Eq. (|) and (||) into Eq. (|^), one gets: 

Zem = \nwc\ ksTln-^^^-^^- ksT^^ . 

(8) 

Because fiwc is negative and much larger than ksT in 
absolute value, one can see from Eq. (^) that, if A is 
not very close to a, the last two terms on the right hand 
side can be neglected. In this case, V(0) and therefore 
a* are positive. Recalling that the bare charge density 
— (7 of the macroion is negative, we see a charge inversion 
effect: the amount of Z-ions adsorbed by the macroion 
surface is larger than needed to neutralize its charge. Eq. 
(^ clearly shows that this effect is driven by correlations 
since the correlation part |/ivFc| of the chemical potential 
is the only positive term in the right hand side. 

It is interesting to mention the limiting case of Eq. 
(^ in which Z-ions are point charges (a = 0). At high 
enough bulk concentration c or low enough temperature 
T, Equation (||) then gives Zeip{0) = \nwc\- For a spher- 
ical macroion, its surface potential is '0(0) = Q* /tmD 
(rj\/ is the macroion radius, Q* is the macroion net 
charge). This givesH: 

Q* = Q.M./QZ'e , (9) 

where Q is the macroion bare charge. This is the max- 
imum possible inverted charge of a spherical ion. Sur- 
prisingly, this charge does not depend on the macroion 



radius. This function Q*{Q) is plotted in Fig. |3|by the 
dashed line. 

For Z-ions with finite radius a, when A approaches 
a the last two terms in Eq. (^) increase substantially 
and reduce charge inversion. When a near complete 
layer of Z-ion is formed [A — a <^ a), these terms 
eventually become bigger than l/iivd- When this hap- 
pens, cr* becomes negative and charge inversion disap- 
pears. If one increases the macroion charge further, the 
second layer of Z-ions appears on the surface of the 
macroion and the charge inversion recovers. Now, one 
can still use Eq. (||) to calculate charge inversion pro- 
vided the bare charge of the macroion is replaced by 
the net charge of the macroion with the first layer of 
Z-ion (Ti = — cr -|- Ze/2^/?>a? . Similarly, as the macroion 
ion charge increases further one obtains a periodic be- 
haviour of charge inversion with respect to the increase 
of macroion charge: charge inversion disappears every 
time a layer of Z-ions is near full and recovers when new 
layer appears. This oscillating behaviour is shown in Fig. 
^ where the solution to Eq. (g) is plotted. In Fig. ||, there 
are small regions near the beginning of each layer where 
the net charge density of the macroion with full layers is 
so small that F < 1. In such a region, Eq. (^) cannot be 
used to describe screening of the macroion surface. In- 
stead, when F < 1 , correlations can be neglected and one 
should use the standard Poisson-Boltzmann mean-field 
screening theory. Hence, in this region the net charge 
is the same as the bare macroion charge (ct* < 0). As 
cr increases, however, F soon reaches unity and charge 
inversion takes over. 

It should also be noted here that, in all above calcu- 
lations, the dielectric constant of the macroion was im- 
plicitly assumed to be equal to the dielectric constant D 
of water solution so that all Z-ions lie on the macroion 
surface. In general, the dielectric constant of a macroion 
is smaller than water so that a Z-ion approaching the 
surface sees an additional repulsive force from its elec- 
trostatic image in the macroion due to the discontinuity 
in the dielectric constant at the macroion- water interface. 
This additional image pushes the Z-ion off at a distance 
from the surface and reduces the attractive force created 
by the correlation hole in the adsorbed Z-ion layer calcu- 
lated above. However, all the oscillating charge inversion 
picture mentioned above remains valid qualitatively. In- 
deed, it is shown in Ref. |ll| that the additional image 
charge lifts each Z-ion off the surface by a distance of 
the order oi A/A and reduces the correlation energy ap- 
proximately by half, thus ^wc remains highly negative 
and charge inversion is preserved. 



III. MONTE-CARLO SIMULATION 



To demonstrate the oscillating behaviour of charge in- 
version, we carry out Monte-Carlo simulations of a small 
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system of a macroion with Z-ions. The Z-ions are mod- 
eled as charged spheres with charge Z = A and radius 
a = 0.9^_B, where the Bjerrum length Ib = e^/DkBT. 
The system is simulated at room temperature (298°K) 
so that Ib — 7.2A. The macroion has radius tm — 3.51b- 
Thus, a full (hexagonal close-packing) first layer of Z-ions 
would contain 47rr|^/2v^a^ ~ 86 Z-ions if one does not 
take into account defects caused by the finite curvature 
of the macroion surface (one cannot put a perfect crystal 
on the surface of a sphere). The actual number of Z- 
ion at full filling should be smaller. In our Monte-Carlo 
simulation, in part also due to the divergence of pressure 
mentioned in the previous section, the number of Z-ions 
in the first layer never exceeds 75. 

The macroion charge — Q is varied from 20 to 550 
so that at maximum there are two layers of Z-ions on 
the macroion surface. Energies of interaction between 
two Z-ions and between each Z-ion and the macroion 
were calculated according to unscreened Coulomb law 
(Ze)'^ / Drij and —ZQe^/DrMi, where is the distances 
between centers of the two Z-ions numbered i and j, and 
TMi is the distances between the centers of the Z-ion 
numbered i and the macroion. The net macroion charge 
Q* = —Q + NZe is determined from the result of the 
simulation where is the number of Z-ions adsorbed on 
the macroion. A Z-ion is considered adsorbed if its cen- 
ter is found within a distance 3a from the top full layer 
of the Z-ion (or from the macroion surface if there's only 
one partially filled layer). The whole system is positioned 
inside a hard spherical shell with large radius L = 17 AIb 
with the macroion fixed at the center of the shell (see 
Fig. §. 




FIG. 5. Schematic drawing of the Monte-Carlo cell. The 
macroion is fixed at the center. Some Z-ions condense on the 
macroion and overscreen it. The rest of the Z-ions spread out 
over the outer boundary of the system. 

In our simulated system, the total charge of the 
macroion and all the Z-ions is 220e, so that the total 
charge of Z-ions is larger than the macroion charge. This 
is needed to observe charge inversion. In practical situa- 
tion, the excess charge of Z-ions is neutralized by mono- 
valent ions at the distance of the order of the screening 
length from the macroion surface. In our simulation, 
these excess Z-ions are neutralized at the external hard 
spherical shell instead (it would be too time consuming to 
simulate the monovalent ions). This difference, however. 



produces no significant impact on the charge inversion 
because of the short range nature of correlations. In- 
deed, it is shown in Ref. |^ that if rg is larger than radius 
of the macroion, charge inversion only weakly depends 
on rg. This is because the capacitance of the macroion 
screened at the distance much larger than its radius only 
weakly depends on r^. In the same way, charge inver- 
sion of the macroion in the shell is almost independent 
on the shell radius when it is much larger then macroion 
radius. This is because the capacitance of macroion-shell 
capacitor is almost constant, while it is charged by the 
negative chemical potential of Z-ion which is determined 
by interactions at the distances of the order of the size 
correlation hole. A, from the macroion surface. 

For most of the Monte-Carlo attempted moves, the 
new position of a randomly chosen Z-ion is chosen with 
uniform distribution in a cube with size Ib centered at the 
old Z-ion position. However, once in every 100 attempts, 
a longer attempted move is made: If the chosen Z-ion is 
found in the vicinity of the macroion {tm < r < tm + 2a, 
where r is the distance from the Z-ion to the macroion 
center), it is attempted to move to the vicinity of the ex- 
ternal shell (L — 3.5a < r < L). Vice versa, if the chosen 
Z-ion is found in the vicinity of the external shell, it is 
attempted to move to the vicinity of the macroion. This 
is done to quickly equilibrate the system and to overcome 
the huge Coulomb barrier which results from the combi- 
nation of short range attraction of a Z-ion to its image 
in the SCL and the repulsion from the inverted macroion 
charge. The usual Metropolis algorithm is used to re- 
ject or accept a move. Typically, in a simulation run, the 
number of attempted moves per Z-ion is 2 x 10^, of which 
the last 1 X 10^ moves are used for averaging oi Q* . 

A snapshot of the complex of the macroion and the 
Z-ions condensed on it is shown in Fig. ^ for the case 
the macroion charge is — 220e. The net macroion charge 
in this case is 12e (the macroion is overscreened) . 




FIG. 6. A snapshot of the macroion with the layer of Z-ions 
condensed on it. The macroion charge is — 220e. AH Z-ions 
found within a distance less than 3a from the macroion are 
shown. In this snapshot, there are 58 of them making the 
macroion net charge 4e x 58 — 220e = 12e. 

The simulation results are shown in Fig. 1^, where the 
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macroion net charge Q* is plotted against the macroion 
bare charge Q together with our theoretical prediction. 
As we expected, the charge inversion disappears {Q* goes 
to zero) near complete filling. After that, the second layer 
of condensed Z-ions starts to form and the macroion is 
undercharged initially. When Q increases further so that 
more Z-ions come to the second layer and correlations 
become important, charge inversion soon recovers and 
increases with Q. 
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FIG. 7. The net macroion charge Q* vs. its bare charge 
Q as the result of Monte-Carlo simulations (the solid circles). 
The solid line is the theoretical prediction using Eq. ( |ll| ) . We 
do not extend theoretical predictions for Q < 50 and for the 
beginning of the second layer because, there, predictions of 
Eq. dllj ) loose their accuracy (see the text). The dash line 
corresponds to the corrected theoretical prediction near the 
full filling, which takes into account lattice defects. 

In order to compare our theory with simulation results 
some modifications are necessary. Because of the pres- 
ence of the hard shell which confines the Z-ions, Eq. (|^) 
should read: 



^ + Zei/'(0) = fcBTln(cLWo) + ZetfjiL) , 



(10) 



where in the right hand side, the bulk concentration c 
is replaced by the concentration cl at distance L. Be- 
cause cl is small, we still assume ideal solution neglect- 
ing Z-ion correlations (so that the chemical potential of 
Z-ions at distance L is kBT\n{cLVQ)). The potential 
ip{L) = ZeQ* / DL obviously does not vanish at distance 
L. 

Correspondingly, Eq. (||) is replaced by: 



Ze[V'(0) - ^{L)] = \^lwc\ - keTln 



1 



A- a 



(11) 



In the simulation, the Z-ion concentration cl (hence, 
their chemical potential) is kept relatively constant by 
choosing the number of Z-ions at any Q in such a way 
that after neutralizing the macroion charge, there are 55 
Z-ions left. As one can see from Fig. 0, the oscillations of 
the number of Z-ions adsorbed at the macroion surface 
do not exceed 4. Such a small deviation in the number 



of free Z-ions from 55 only weakly affects c^. Further- 
more, due to the excess charge of these extra Z-ions, they 
spread out in a thin layer near the outer boundary where 
their charges are neutralized. This helps to keep cl rel- 
atively constant. Because in Eq. (^), cl appears under 
the logarithmic function, the chemical potential of Z-ions 
oscillates even weaker than cl- 

The solution to Eq. ( [TTI ) is plotted by the solid line in 
Fig. 1^. As we can see, although the theoretical prediction 
agrees favourable with the Monte Carlo results, it overes- 
timates the net inverted charge near the full filling. This 
discrepancy is due to the use of a planar hard disc system 
in the theory of the Z-ion layer at the macroion surface. 
Indeed, unavoidable lattice defects caused by the curva- 
ture of the macroion surface (there are 12 disclittations 
and a number of dislocations in the Z-ion layerES) can 
lead to a large misprediction of the number of Z-ions at 
complete filling if we assume a flat layer with the same 
area. For example, calculations using icosahedron show 
that even in the optimal situations with "magic" num- 
bers of Z-ions, planar geometry mispredicts the number 
of Z-ion at complete filling by 1 for iV = 12 and by 3 for 

= 32. Away from the magic numbers where a number 
of dislocations appear, this misprediction is even larger. 
The dash line in Fig. ^ is the adjusted theoretical pre- 
diction of Eq. ([HI ) assuming the number of Z-ions at 
maximum filling of the first layer is 79 instead 86. The 
agreement between this curve and the Monte Carlo re- 
sults is much better than the solid curve. It should be 
noted that, defects are important only near complete fill- 
ing where the Z-ions layer forms a real solid. At smaller 
filling where Z-ion layer is liquid like, defects play no 
role. That is why the adjusted prediction (the dash fine 
in Fig. ^ is not drawn at smaller filling. Obviously, at 
larger macroion size, the region where defects are impor- 
tant is relatively smaller. 

The planar approximation of the macroion surface 
also forces us not to extend theoretical prediction to 
Q < 48 (or 12 Z-ions) because, at small Q, the av- 
erage distance between Z-ions on the macroion surface 



2A ~ y8'!T{rM +a)^Ze/V3Q is larger than the radius 
rj\/ of the macroion. This makes the concept of a two- 
dimensional Z-ion layer, and therefore Eq. (|l l|) , quanti- 
tatively incorrect. Again, if the macroion size gets larger 
or the Z-ion size gets smaller, the Z-ion layer can be 
truthfully regarded as two-dimensional in a larger region 
of Q. This in turn makes Eq. ( [Tl| ) applicable quantita- 
tively in a wider region of Q. It should be noted that, 
even though Eq. (^l|) fails at small Q in our case, the 
qualitative picture of charge inversion is still valid. This 
is because, in our case, the plasma parameter F is still 
large and does not reach unity until (5 — 2. 
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IV. VOLTAGE CONTROLLED OSCILLATIONS 
OF CHARGE AND CAPACITANCE 



In the previous sections, we discussed the osciUating 
behaviour of the inverted charge as a function of the 
macroion charge. GeneraUy speaking tuning the charge 
of an insulatiag macroion is not a trivial task, although 
can be donen3. We propose here an alternative exper- 
iment that would show this oscillating behaviour when 
the macroion surface potential is tuned instead. It is 
sketched in Fig. W. 
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FIG. 8. A sketch of the proposed experiment. Two metal- 
lic plates are separated by an insulator (white) and a water 
solution (gray) containing multivalent Z-ions (shown by solid 
circles). Under the voltage V applied to the capacitor, charges 
— Q and Q are induced in the metallic plates. Z-ions with to- 
tal charge Qt = Q* + Q are adsorbed on the right side of the 
insulator. The excess charge Q* due to charge inversion is 
linearly screened at the distance rs by monovalent salt. 

The role of the macroion is played by the left metallic 
plate covered by an insulator layer with thickness d, sur- 
face area S and dielectric constant Di. The other side of 
the dielectric is in contact with the ionic solution which 
contains the Z-ions. Some of Z-ions condense on the di- 
electric surface to screen the electric field created by the 
negative charges on the metal plate. By changing the 
voltage V applied to the other plate, one can easily vary 
the charge of condensed Z-ions on the dielectric surface. 
We show below that if one measures the differential ca- 
pacitance of this systems as a function of V , one can 
expect an oscillating behaviour. 

To deal with a macroion at a constant potential, our 
theory needs some modifications. When the voltage V 
is applied, a charge —Q is induced on the left metallic 
plate. Z-ions layer with total charge Qt is then formed 
on the dielectric surface in solution to screen the electric 
field of the plate. The excess charge Q* = Qt — Q due 
to overscreening is in turn linearly screened at the dis- 
tance in the solution. Due to the electro-neutrality of 
the solution, charge —Q appears on the other side of the 
solution container to screen the second metallic plate. If 
d is large enough {d ^ A) so that one can neglect the 
image of Z-ions in the left metallic plate, the free energy 
of the system can be written as: 



2C(d) 2C(r,) 



+ Ne{n) + 



-QV- Nfib , 
(12) 



where Q is the charge, C{d) = DiS/And is the capaci- 
tance and Q^/2C{d) is the energy of the planar capac- 
itor formed by the dielectric layer with thickness d. As 
we mentioned above, the excess charge Q* is linearly 
screened at distance in solution. The free energy of 
this screening atmosphere can be expressed as the energy 
of a capacitor with charge Q*, thickness and capaci- 
tance C{rs) = DS/AttTs (The second term in Eq. (p^). 
The third term in Eq. (n2|) accounts for the correlation 
energy of the SCL of Z-ions condensed at the dielectric 
surface [N = Qt/Ze is the number of Z-ions condensed 
and n = N/S is their two-dimensional concentration). 
The fourth term in Eq. ( p^ ) is the energy of the capac- 
itor with thickness on the right side of the solution 
container. On this side only monovalent ions screen the 
electrode. Therefore correlations are weak and screening 
can be described by the usual linear theory which gives 
the free energy Qy2C{rs). The fifth term in Eq. (|l|), 
—QV, is the work the source does to maintain the system 
at the constant voltage V. Finally, the last term in Eq. 
( p^ ) is the change in the free energy of the bulk solution 
which acts as a reservoir to provide the Z-ions to the 
macroion surface. 

Minimizing the free energy ( p2| ) with respect to the 
unknown parameters Q and Qt, one gets: 

V^Q {C(d)-' + 2C(r,)-i) - QtC(r,)-i ^ (^3) 
-(mH - ^ib)/Ze =. Q*C(r,)-i . (14) 

As shown in previous sections, due to correlations be- 
tween Z-ions, —{ii{n) — fib) is positive everywhere except 
near complete filling. According to Eq. (14), this means 
that Q* is positive, indicating charge inversion. 
Solving Eqs. (|l|) and one gets 



V = Q{Cid)-' +Cirs)-') + 



fi{n) - fib 
Ze 



(15) 



which shows that the correlation between Z-ions con- 
densed on the dielectric surface works as an additional 
voltage applied to the system. 

A measurable quantity of the system is its differential 
capacitance C — dQ/dV, or 



c- 



dV 
dQ 



C(d)-i -f C(r,)-i + 



1 dfi{n) 
Ze dn 



(16) 



As we saw in previous sections, excluded volume effect 
leads to the oscillations in the chemical potential fi{n) of 
the Z-ions which becomes positive every time a full layer 
of Z-ions is developed. From Eq. (p^, this ultimately 
results in the oscillations of the inverse capacitance C"^ 
of the system. Thus measuring C{V) one can indirectly 
observe the oscillations in charge inversion. 
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Concluding this section, we would like to mention that 
the oscillations of the differential capacitance we obtained 
above are remarkably similar to the well known oscil- 
lations in the compressibility and magneto-capacitance 
in thej-two-dimensional electron gas in quantum Hall 
regimec3. In this case, the oscillations are related to the 
consecutive filling of Landau levels. The compressibil- 
ity of the two-dimensional electron gas vanishes every 
time a Landau level is fully occupied. In our system, 
the two-dimensional electron gas is replaced by the quasi 
two-dimensional liquid of Z-ions and the magnetic field 
is replaced by the varying voltage. In this sense, we can 
say that we are dealing with a classical analog of the 
quantum Hall effect. 



V. CONCLUSION 



creases until, at some value of a which the authors call 
(Te, the blobs start overlapping with each other. In con- 
trast with the hard-core repulsion case, blobs are easily 
compressible. Therefore, when the surface charge a in- 
creases further, the adsorbed PE molecules form a scaling 
structure of layers of blobs with increasing size and de- 
creasing density from the bottom to the top. In the top 
layer, blobs have the same size as unperturbed blobs in 
the solution. As in the case of hard core compact .Z-ions, 
it is this top layer which is responsible for the image of 
an approaching new PE molecule and, therefore, for the 
correlation induced charge inversion. According to Ref. 
p!^ , in the case of flexible weakly charged PE the proper- 
ties of the top layer do not oscillate with a because newly 
adsorbed PE molecules lead to an additional compression 
of lower layers. Thus at cr > cTg, the inverted charge does 
not oscillate and remains the same as at cr = de. 



Summarizing our results, we have shown that the ex- 
cluded volume of .Z-ions can dramatically change the 
screening of a macroion by Z-ions. In particular, charge 
inversion is shown to oscillate with the macroion sur- 
face charge density a. At small macroion surface charge 
density, it increases with a. However, when the surface 
charge is high enough so that a near full layer of Z-ion 
is formed, charge inversion disappears quickly because 
of the vanishing compressibility of a full layer of Z-ions. 
Charge inversion recovers when a increases further so 
that a second layer of Z-ions builds up. Since tuning 
of the charge of an insulating macroion is not a simple 
task, we also propose an experimental setup with metallic 
electrodes playing the role of a macroion. This experi- 
ment can be used to indirectly observe these oscillations 
of charge inversion. 

Oscillations of charge inversion discussed in this paper 
are related to the hard-core repulsion of spherical Z-ions. 
Similar effect should exist for the adsorption of rigid rod- 
like polyelectrolytes (PE) on oppositely charged surface. 
On a strongly charged surface rods can form more than 
one layer leading to a deep minimum of the charge inver- 
sion every time when the top layer is completely full. It 
is important that charge inversion in both cases is deter- 
mined by the top layer. 

It is interesting to discuss the relation between our 
work and the recent theory of adsorption of a thick layer 
of weakly charged, flexible PE with excluded volumetil in 
which no oscillations of charge inversion with the increas- 
ing surface charge density of macroion, a, was found. 
Instead, the inverted charge reaches a constant maxi- 
mum value when PE molecules start to overlap. The 
picture described in Ref. |lj is as follows. A flexible PE 
molecule in solution can be viewed as a rod-like string 
of Gaussian electrostatic blobs. Each such blob has a 
quite low density of polymer. At small absolute value 
of CT, PE rods screen the macroion surface forming cor- 
related liquid of locally parallel rods which overcharges 
the macroion. As a increases, charge inversion also in- 
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VI. APPENDIX 

To derive the chemical potential (^, let us start with 
the standard relationship between the free energy and 
the canonical partition function. 

F(iV,t/,T) = -fcsTlnZ (17) 

where the canonical partition function Z is given by: 



Z = 



1 



,N 



N\ 



dri...drN exp < — 



1 



■ N 

E 

.4=1 



Zeve{ri)+ 



N 

E 



[vc{\ri - Tjl) + Vhci\ri 



(18) 



Here Ve{ri) is the electrostatic external potential at the 
position of the particle i (in our case, Ve{r) is the po- 
tential due to the macroion charge), vc{r) = Z^e^ / Dr is 
the usual Coulomb part of the pair potential and Vhc{r) 
is the hard-core repulsion part: 



Vhc[r) 



oo if r < 2a 
if r > 2a 



(19) 



and So is the normalizing area. In writing down Eq. (|l^) , 
we use the primitive model of the solution where all wa- 
ter degrees of freedom are averaged out, resulting in an 
effective dielectric constant D and an effective quantum 
cell area sq. 

Substituting Eq. (^ into ( p^ ) and differentiating with 
respect to the particle charge e, one gets 
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dF 



AT 



TV 



i<3 = l 



(20) 



here the notation (...) denotes the statistical average. 

From Eq. (pO|), one gets the final expression for the 
free energy of the system: 



F — Fq 



2de 



N 



N 



\i=l i<j=l I 



(21) 



where Fq is the free energy of the systems of hard discs 
with radius a at the same concentration. The integration 
denotes the electrostatic contribution to the free energy. 
As in the case of a SCL of point charges, because V is 
large, one can approximate this electrostatic contribution 
by that of a Wigner crystal: 



F ~F^ + Newcin) 



(22) 



Due to hard-core repulsion, the Z-ions are never at a 
distance smaller than 2a from each other. Therefore, the 
approximation made in Eq. ( |2^ ) is even better than that 
of point charge SCL and is asymptotically exact when 
A—> a. 

Eq. (M) can now be easily obtained by differentiating 
Eq. (E2j with respect to the number of Z-ions N . 
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